Local semicircle law in the bulk for Gaussian ^^-ensemble 



Philippe Sosoe Percy Wong 



Abstract 



We use the tridiagonal matrix representation to derive a local semicircle law for Gaussian beta 
ensembles at the optimal level of n~^'^^ for any 6 > 0. Using a resolvent expansion, we first derive a 
semicircle law at the intermediate level of n"^/^"*"*; then an induction argument allows us to reach the 
optimal level. This result was obtained in a different setting, using different methods, by Bourgade, 
Erdos, and Yau in 3 and in Bao and Su _2_. Our approach is new and could be extended to other 
tridiagonal models. 



1 Introduction 



Given a symmetric n x n matrix M, the Stieltjes transform sj\/ of its eigenvalue distribution function, fj,M, 
is defined by: 

Ss{z) = / — — /iM(dx), 



Mm((— oo,a;]) = — • tJ{A < a; : A is an eigenvalue of M}. 
n 

The Stieltjes transform is an important tool in the study of asymptotic properties of random matrices. 
Several proofs of Wigner's theorem on the eigenvalue distribution of random symmetric matrices in the large 
n limit are based on identifying a pointwise limit of the Stieltjes transforms with the Stieltjes transform of 
the semicircle density: 

2 , 

Psc{dx) = - VI - dx. 



Up to a constant, the imaginary part of the Stieltjes transform is the convolution of fiM with the Poisson 
kernel: 

3ss(z) = nP,f * iis{z) =TT — - HM{dx), E + if]. 

J [X — hy + 7]-^ 
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Hence the imaginary part Ssm is a smoothed version of the eigenvalue density, and the imaginary part 
r] = Q'z appears as a parameter controlhng the resolution of the approximation. In their work on universality 
for Wigner matrices (see [S], [TJ), Erdos et al. developed various versions of the local semicircle law. Let 
Sn be the Stieltjes transform of a normalized Wigner matrix W: 



Then s„(z) and Ssciz), the Stieltjes transform of the semicircle distribution, are asymptotically arbitrarily 
close as long as ^ For precise statements, see [TU], [TT]. The local semicircle law is a crucial input 
in proofs of universality for matrices with independent entries. Information on s ioi z = E + iri translates to 
information on the empirical distribution on scale 77, i.e. control of distribution of the eigenvalue counting 
functions of intervals of length of order 77' > 77: 



Away from the edges of the support of the semicircle density, such intervals typically contain nrj' eigenvalues. 
See for example Lemma 16.31 below. 

Apart from Wigner matrices, the central objects of interest in random matrix theory are unitarily invariant 
matrix ensembles, with joint eigenvalue density of the form: 



with /? = 1, 2, or 4. For Gaussian potential V{x) = (/3/4) • x^, this is the eigenvalue density for GOE 
(/? = 1), GUE {f3 — 2) and GSE (Gaussian Symplectic Ensemble, (3 = 4) matrices. For general /? > and 
V, the formula above defines the /3 ensemble density. In [3], the authors compare the local statistics of n 
"eigenvalues" distributed according to the beta ensemble density with a convex, analytic potential V{x) 
and general /3 > to the statistics of the Gaussian ensembles. They show that the statistics of /3 ensembles 
are universal in the sense that they asymptotically match those of the Gaussian ensemble with the same 
value of /3. Although the proof of universality for (3 ensembles in [3 differs from those in the Wigner case 
in many respects, the initial step is also to derive a fine localization result for "eigenvalues" of /3 ensembles 
with convex analytic potentials. The convexity assumption was subsequently removed in 0]. 

The purpose of this paper is to give an alternative derivation of a local semicircle law for Gaussian (3- 
ensembles. Our approach is based on the tridiagonal matrix model introduced by Dumitriu and Edelman 
[8] , and thus essentially different from that in 3 . Note that Bao and Su [2j have derived a local semicircle 
law, down to scale l/^/Togn, also based on the tridiagonal representation. Their method does not provide 
an estimate on the probability of deviation other than that this probability vanishes as 7^ — >■ c». 

The main result is the following: 

Theorem 1.1. Let Ssc be the Stieltjes transform of the semicircle law, and s be the Stieltjes transform of 
the normalized Gaussian (3-ensemble measure (V — jx"^ ). Then, for any c, fc, e > and < 6 < 1, there is 





const. X 
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a constant Cc,k,e such for z = E + irj in the region 

D^,s ■■= {z:^z> r^-l+^ -I + 5 <^z < I - 5,\z\ < 10}, 

we have: 

f{\ssc{z) - SnAz)\ > c) < CcKc^sn-''. 



The following corollary can be derived from the theorem, see e.g. [3], [TU] 
Corollary 1.2. Let the semiclassical location be defined as the number such that 

[ Pscis)ds = ^ 

then for any 5,e > Q, n > 1 and any j £ [6n, (1 — 5)n], 

|A,-7,| <n-i+^ 
with probability no smaller than 1 — Ck.e.S'n^'' ■ 

Remark 1. The availability of a tridiagonal matrix model is not specific to the Gaussian case and our ap- 
proach can potentially be extended to give local semicircle laws for other matrix ensembles with tridiagonal 
representations, such as the Laguerre ensembles considered in [8]. 



The entries of the tridiagonal representation of /3 ensemble matrices have very different sizes. The approach 
based on Schur complementation and predecessor comparison used in [lOj . as well as further works on 
the local semicircle law relies on the self-similar structure of Wigner matrices and thus cannot be directly 
applied in our setting. 

As explained in Sectional we proceed by separating the model matrix into its expectation, a deterministic 
Jacobi matrix whose resolvent can be analyzed using asymptotics for orthogonal polynomials in the complex 
plane (Section [S]) , and a fluctuation part, whose contribution is shown to be asymptotically negligible. To 
control the resolvent expansion, we require some estimates for the entries of the resolvent of the expectation 
matrix; these are established in Section [S] The expansion provides us with control of the resolvent down 
to 5z > n~^/'^~^'^: an inductive argument in Section [5] establishes the local semicircle law in the region 
3z > 



2 Notation 



Definition 1. Let /3 > 1, the Gaussian /3-ensemble is an ensemble of eigenvalues in M" that have the 
following probability density: 

d/(Ai,...,A„) = G„,^ H |A,-A,f exp[-^f^An[]dA„ (1) 

l<i<j<n \ 1 / i 
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where Gn.p is a normalization constant. 

Definition 2. Given a measure ^ supported on R, the Stietjes transform of ^ is defined to be the complex- 
valued function 

s^{z)^ [ ^—dfiix). (2) 

J X - z 

The Stietjes transform is analytic on the upper half plane and converges weakly to fJ.{x) as 5z — >■ . 

Definition 3. The Hermite polynomials Hn{x) are orthogonal polynomials with respect to the weight 
e~^^ on the real line (this is the "physicist's normalization"). The n-th Hermite polynomial has leading 
coefficient 2" , and is given by the Rodrigues formula: 



The n-th Hermite function En is defined as: 



The functions En, n > form an orthonormal set in L^(M). 



Dumitriu and Edelman 8 , generalizing earlier observations by Trotter [18 in the case {3 — 2, introduced 
the tridiagonal model for the Gaussian /3-ensembles, among other matrix models: 



Theorem 2.1. Consider the matrix given by 



An,l3 = 







/ iV(0,2) X(n- 

X{n-l)l3 ^(0,2) X(n-2)/3 

X(n~2)P N{0,2) X(n-3)/3 



V 



X3/3 





^(0,2) X2P 
X2f5 N{0,2) 
Xp 









X/3 

N{0,2) J 



(3) 



where N{0, 2) denotes a random variable whose distribution follows the Gaussian distribution with mean 
and variance 2; Xk is a random variable having a chi distribution with k degrees of freedom. The upper 
triangular part of the matrix consists of independent random variables and the matrix is symmetric. The 
joint density of the eigenvalues of An,p coincides with the fi-ensemble density. 



We will refer to An,p as the symmetric model for the /3-ensemble. Conjugation by the diagonal matrix 

Dn,p = diag(di, . . . , d„), (4) 
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where 

i-l 

d,(n,/3) = (/3n)(i-^)/2[]x(„-j)/3. 

(see e.g. shows that df is also the joint probabihty density for eigenvalues of the (non-symmetric) 

matrix 



A. 



N{0,2) 






A^(0,2) 







/]3^ 







iV(0,2) 




7V(0,2) 



X0 









iV(0,2) / 



Throughout the rest of this paper, we will be switching between the symmetric and the asymmetric models 
and will keep track of the changes involved. 



For the rest of the paper, we shall use the following notion of an event depending on some index n having 
overwhelming probability: 

Definition 4. We say a sequence of events En holds with overwhelming probability if for all n, P(-E„) > 
1 — Oc{n''^) for every constant C. 



It should be observed that a union of events of overwhelming probability for some fixed fc > still holds 
with overwhelming probability. 



Lastly, C > will always denote a constant whose exact value need not concern us. 
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3 Local semicircle law for zero temperature 



Let 71 be a positive integer, and consider the "zero-temperature" /3 ensemble matrix: 



An 



1 

71 



/ VrT^ 





V 










V2 










1 





1 


) 



(6) 



This corresponds to the /3 — )■ oo hmit of the symmetric model matrices An^p ■ 



We denote the trace of the resolvent Soo ■ C+ C+ of An. oo by 

, , 1 / 1 . 1 ^ 1 

Soo(Z) — — tr —=An oo ~ Z \ — ~ y , T ■ 

n \ J n \j — z 

The rescaled resolvent converges uniformly on compact subsets of z \ C (see [5], p. 159-167): 

Here Ssc is the Stieltjes transform of the semicircle density: 

Ssc(2) = - / dx. 

TT a; — z 

As a first step towards the proof of the local semicircle law for /3 ensembles, we establish the following 
quantitative version of this result, close to the real axis: 

Lemma 3.1. Let < (5 < 1. There exists a constant Cs,e such that, for any z with 1/10 > 3z > n~^^'^ , 
and —l + 5< 5Rz < \ — 6, we have: 

\sAz)-Ssc{z)\<Cs,,n-\ (7) 



The proof makes use of the "Plancherel-Rotach" asymptotics for Hermite polynomials in the complex plane 
obtained by Deift et al. [5] using the Riemann-Hilbert approach. Their methods yield uniform error bounds 
which imply our result. 
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3.1 Proof of Lemma 13.11 



The eigenvalues Aj, 1 < j < n of A„ oo are given by the n distinct, real zeros of the n-th Hermite polynomial 
Hn, rescaled by ^/2n, of. [7], [S]. For any z, the resolvent s{z) can be written as 



1 " 1 



-1 ^'3 



A, — z 



Hn{V2^z) HniV^z) ' 

where the prime denotes differentiation. The derivatives of Hermite polynomials satisfy the identity: 

H'Jz) = 2nHn-i{z). 
Hence, the trace of the normalized resolvent of A„^oo has the expression: 

We will check that the final expression is close to the Stieltjes transform of the semicircle density. 

In [HI , the authors derive asymptotics for general orthogonal polynomials with respect to even polynomial 
weights by means of a rescaled Riemann-Hilbert problem. For z as in the statement of Lemma l3.1[ formula 
(8.32) and Theorem 7.10, and the remarks in Appendix B. in [6 imply: 



g„(^/2^z) ^_„^(,) ^ 1 f {z-lf/^ ^ {z + lf 



1/4 



23n/2„n/2 2 \{z + lY/^ (z- 1)1/4 



1)1/4 (^_,_ 1)1/4 



2^ \{z + lY/^ (z- 1)1/4 



0(l/n), 



with an error term uniform in the region of interest. The function g is the logarithmic potential of the 
equilibrium measure associated to the Hermite polynomial: 



2 /"i 

-jiz) — — I \/\ — logfz — x) da;, 
71" J-i 



defined for z e C \ (—00, 1]. The function is given by: 

(p(z) = ~iz{\ - z)i/^(l + z)i/^ - 2% arcsinz + itt. 
Here the square roots are principal branches, e.g.: 

(1 - z)i/2 = exp((l/2) log |1 - z| + i(l/2) arg(l + z)) 
where arg e (— tt, tt). The function arcsin is the inverse of the conformal mapping 

sinw : \\^w\ < ^| ^ C \ (-cx),-l] U [l,oo). 
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Note that smw maps {Qw > 0} one-to-one and onto the upper half-plane {5z > 0}. 

The Riemann-Hilbert problem also provides us with asymptotics for the rescaled n — 1st polynomial. 
Indeed, computing the second row of the solution to the Riemann-Hilbert problem and using Stirling's 
approximation, one finds (cf. (3.9), (5.54) and (8.32) in 'B^): 



2(3n-l)/2„(«-l)/2 2V (Z + 1)1/4 (2-1)1/4 



1/4 



^ > ' ' > ' Oil/n). 



(Z + 1)1/4 (z^ 1)1/4 

The error term is 0(l/n) by the uniform boundedness of (z + 1)^^/^ and (z — l)^^/^ for the values of z 
that concern us. 

We now argue that the factor e""'^'^-' is rapidly vanishing in n. We have: 

^ip(z) = sRz • |1 + z|i/2|i _ z\^/^sm ((l/2)(arg(l + z) + arg(l - z))) 

-I- 5z • |1 -I- z|i/2|i _ z|i/2cos ((l/2)(arg(l + z) + arg(l - z))) + 23arcsinz. (9) 

The last two terms are always positive. For 3z > 0, — 1 + (5 < 5Rz < 1 — (5: 

arg(l-z) e (-7r/2,0), 
arg(l + z) e (0,7r/2). 

Note also that > for 5Rz < 0. We will need a lower bound for the real part in the region where is 
positive. For such z, 

\Z - 1|1/2|1 + z|l/2 < (1 _ g?z)l/2(l -f 3?2)l/2 (1 + (9z)2)^/' . 

On the other hand, since the argument of the sine is negative for < ^z < 1 — ^: 

1 / 5z \ 1 / 3z 

sm(arg(l + z) + arg(l ~ z)) > - arctan I — I — - arctan 



2 \l + ^zj 2 \l-^z 

> 



(l + 5Rz)(l-3fiz)' 



We compare this to a lower bound for Qw = Sarcsinz. By our definition of arcsin: 



Qz = COS^W, (10) 

3?z = sinSftw. (11) 
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For 3fJz > 0, this implies cos^w, sinSRw > 0. Using (ITU)) . ([TT|) . wc find: 

sinSRui < 5R2 

cos^w > (1 - + Sff^)l/2 

Using (|10p again, this imphes 

5u; < (2<5)-i/25z. 

Without loss of generality, we may assume 5z < 5/100. Now using = 1 + e'^x for < a; < c, we have: 

2 < e'-*"" + e"^"" < 2 + e^/i^Gu;, 
3?z • (1 - (ei/i°/2) • ^w) < sinSRw, 
29t(; < e'-**" - e"^™ < 23w + 

cos3?ix; < (1 - i^zf))^^^ + lOQw ■ \2^z - e^^/^°Qw\. 

Hence we have: 

(1 - 5Rz)i/2(l + ^z^^^^w > 3z + 0(3w)3/2. 
Here O stands for some terms of higher order in multiplying small constants depending on S. 

Putting all the above together, we find: 

A2\ . c> 



_ (2 - {^zY) ■ ^z 

Recalling that 3z < (5/100, we find 



uniformly in 3z > n 
The factors 



(Z- 1)1/4 ^ (^ + 1)1/4 



(0 + 1)1/4 (2 _ 1)1/4 

are bounded uniformly in the region specified in the theorem, so the above gives a sub-exponential decay 
rate for the factors multiplying e""'^'-^-'. 

Using the approximations above, a calculation shows that, uniformly in 1/10 > '^sz > n~^~^'^: 

S„(Z) = 2 • (-Z + (Z - 1)1/2(Z + 1)1/2) ^ 0(1/7^). 

We once again take the principal determinations of the square roots. It is readily verified that the boundary 
values of the imaginary part of the first term on the right for 3z — > 0"*" are given by n times the semicircle 
density. Since at infinity we have 

(Z- l)l/2(z + 1)1/2 + 0(1), 

for our choice of the square roots, it follows that the analytic function defined by the expression on the 
right above is equal to Ssc- 
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4 Semicircle law for Gaussian /3-ensemble at level n^/^ 



In this section, our goal is to prove the fohowing semicircle law at a suboptimal level of n 1/^+*^ for any 
e > 0: 

Proposition 4.1. Let s„(z) be the Stieltjes transform of the measure induced by the eigenvalues of the 
normalized Gaussian (3-ensemble, -^An^p. Let Ssc{z) be the Stieltjes transform of the semicircle law. For 
any r > 0, there exists a constant C{t,j3) independent of n such that, with overwhelming probability, we 
have: 

sup |s(^)-s.c(^)|<C(e,/3)n-^/ioo, (12) 

where for e > and < (5 < 1 the domain D is defined as 

D^^s {z : d > > n^^/^+e^ -1 + d < < 1 ~ S}. 



The idea of the proof is to expand the Green's function around the zero temperature case and estimate 
the differences between the two. Using the asymmetric tridiagonal model ([5]), we can write 

-^An,l3 = ^i„.oo + A (13) 

where A = A(n, /3) is a bidiagonal matrix whose entries are independent iV(0, variables on the main 

diagonal and independent random variables with distributions 2^xfn-k)p ~ ^(2n^^(Ti-fc)/3) main 
sub-diagonal. 



Using the resolvent expansion, we can write 

m 

R\z) = k^{z) + ^(i?°°(z)(-A))?'i?°°(z) + (i?°°(z)(-A))™+ii?'3(z). (14) 

Here, BP(z) — (-i=A„^^ — zL)^^ is the resolvent matrix for the rescaled asymmetric model (the superscript 
is not to be confused with taking powers) and we suppress the subscript n when no confusion is likely. We 
shall also suppress z in future equations when the dependence on z is understood. Taking traces on both 
sides of the previous equation and normalizing by , we have 

sfi{z) = Soo(^) + l^fj tr((i?°°(-A))^'i?°°) j + n'^ tr ((/Z°°(-A))™+ii?'^) . (15) 

The proof of Proposition 14.11 depends on the following estimates on the elements of the resolvent of the 
symmetric matrix ©, 

Proposition 4.2. Let e' = e/8. For each 1 < k < n and 3z > n^^^'^'^'^, we have the estimate 

<min{Cr^l/2(logn)(l + |^/fc- V^gfJz|)-\ni/2-^}. (ig) 
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For k such that \^/k — Ey/ri\ < with < rj < 1/2, we have 

\RTk\ < CE,,n^/'-'' . (17) 

For k ^ I and > n^^/^+e 

\Rm\ < n^^^~'' k-^/H-^/^\Vk - Vl\-\ (18) 



We will first show how to prove ProDOsition l4.1l froni Proposition l4.2l and then proceed to prove Proposition 
I in the next section. 



Proof of Proposition \4-l\ Let us first establish a simple lemma 
Lemma 4.3. 

f^|i?^,r <Cn^-™^' (19) 
fe=i 

for m > 3. 

Proof of lemma. We split the smn according to the distance h — /i(5ftz) — |\/fc — y/n^z\: 

By (|17p . the first term in (PH)) is bounded by: 
To bound the second term, we use (fTSl) to find: 



Cn^(logri)™ E (l + |Vfc- V^SRzl)"™ < Cn'5-(log7i)™ni/-n 



fc:/i>ni/4 



n 

Il,J2,---,Jp 



□ 



A term in the resolvent expansion (jl4p is of the form (suppressing the superscript oo) 



where i'^ takes nonzero value only for i'^ = is or i'^ — is — 1. We want to rewrite this sum in terms of the 
resolvent R of the symmetric model. The asymmetric model matrix is obtained from the symmetric model 
by the conjugation: 
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The transformation rules for the resolvents follow from this; on the diagonal we have 

Rkk = Rkk 

for every k, and on the off diagonal: 

11^ = 1 X(n-j)fj 

As a result, when i'^, = ik in some summand in (j2ip . the products in the numerator and denominator of 
the successive factors and Rii^i^^^ cancel each other out. If i'j. — — 1, after cancellation, there 

remains a factor of the form: 

Xn-k-l 

y/np 

Note that is at most of order 0(1), with Gaussian tails. Thus we may rewrite (pij) as 

i ^ R,,,,rn^'^A,,,,R,,,,fr^'^\,,^...X^,^R,^,,. (22) 

ii,i2,...,ip 



The offdiagonal entries of A are centered x^-square variables with variance 

VarA.,._,.!i^^.O(l/.). 

Since A^; has exponential tails for all fc, with overwhelming probability, we have for any < c < e, 

maxlAfejl < n-i/2+c/4 

k,l 

maxmfe < n ' 

k 



Therefore with overwhelming probability, the sum ([22]) is bounded by 

ii,...,ip 

We are going to use the estimates from Proposition 14.21 The estimates for Rki and Rk-i.i only differ by 
some constant. When k,k — 1 I, this follows immediately from expressions (fT6|) and ([T8|) . When k — I, 
notice that the right side of (|18p is of order , and we certainly have \Ru\ < Cn^^^~'^. Conversely, 

writing Rn-i as a sum over eigenvectors as in the proof of Proposition 14.21 in the next section, and using 
Cauchy-Schwarz and (fTH]). \Ru-i \ can be bounded up to a constant factor by the right side of (fTB]). Thus, 
the above sum is bounded with overwhelming probability by 

Cf„-in'(p-i)/2+(p-iW2 ^ \R,^,J\R,^,J...\R,^,^\, (23) 

for some C > 0. To estimate the sum, first notice that we have the following estimate, holding uniformly 
in j: 

^Vfc^'lTfc- V7r' < (^logn. (24) 
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In the sum, whenever z; 7^ for some I, we gain a power of n~'^ using the estimate and from (|24p . 
the sum over such pairs introduces a factor bounded by Clogn. Whenever ii = and the power m of 
I™ is at least 3, we use lemma (|4.3|) to gain a power of . Due to the repeated index, the power of 
ii in the sum is — 1 and so the sum over ii will be bounded by log n as well. 



To make the last paragraph precise, consider the pairs of indices (ii, 12), («2, is) ■ • • , {ip, *i) as edges on a 
graph (of n vertices). We call an edge exploratory (E) if ii ^ and stationary (S) otherwise. Each 
p-tuple of pairs ("path") belongs to a category (Ci, C2, . . . , Cp), where each C; G {E, S} denotes the type 
of edge (i/, i;+i). A category of the form 

{ E, E,^. .,E^ , S, S,y.,S^ , E,..^.,E^ , • • • , E, ■ - ^'"J ' ^- ^ 

ji times ki times J2 times Jk times times 

corresponds to a partial sum of of the form: 

(25) 

The summation is over the set {1 < < n,ij ^ ik for each j 7^ k}. We will bound the 

contribution to the sum (|23p from each category using (|4.2I) . We will sum successively over all indices, 
beginning with ii. Collecting all factors depending on ii, we find that the sum over this index is: 

^ ^ 1^*1*2 I + + I I I 

»l/*jl + ...jjc + l.'2 

where fc/f > 0. Inserting the estimates and (|18p . we that the last sum is bounded by: 

X V f-i/^ ^ / 1 (26) 



Note that •\/*ji+...jjf +1 1 ^ + ^ jx+i' ^^^^ bound, and perform the sum over ii using 

when kji < 2, we obtain: 

CV/2-^''.(logn)2.z-i/^.z^Yt.,,+i. 

When fc/f > 3, performing a dyadic decomposition around E^/n as in the proof of Lemma [4.3) and using 
on each dyadic piece, we obtain the bound 

CV/-(^-^'=-)^'.(logn)^.z-^/^.^/t..+,.+i- 
We then sum successively over 12,13, ■ ■ ■ ,iji+...+jK ■ For each index, we encounter a sum of the form: 



13 



where k >0. Using (|16l) . ([TS]). and ([M]) as above, this is bounded by: 
if < fc < 2, and 

C'^+ii-Y'"^''''^^'""^'^'^'' log"' 

whenever fc > 3. When we reach the final index after cancelhng the factor Jk+i carried 

over from the summation over ii, we have to sum 

E|E>. \kK-l 

'31 + --.+3jf + 1 

When fcif_i = 0, this is of size n. When kx-i > 3, it is bounded by n^/'^~''>'K-i'^ . 

In summary, for any block of three consecutive edges in any category, we gain a power of — e'. Multiplying 
the contributions from the summations over the p indices, we find that the expression in (|25p is bounded 

by 

C2p^-(p+l)/2+(p-l)c/2 . ^(p+l)/2^-pe73_ 

The last expression is 0{n~P'^^^^) provided c is chosen sufficiently small. Since there are at most 2^ 
categories contributing to (|23|) . all terms in the first sum in the expansion (fT5|) vanish as n — oo. Taking 
m in (|14p large enough, depending on e, to make pe/48) < — 1/2 when p > m, we can bound the remainder 
term of the resolvent expansion using the trivial bound for R^-: |i?fj < n^^^"*^. 



Returning to the sum in the expansion (jisp . the above implies that 

The proposition then follows from the above estimate and Lemma 13.11 □ 



5 Proof of estimates for R 



ki 



In this section, we prove the estimates for i?^ in Proposition reproduced here for the reader's conve- 
nience: 

For $52 > n"^/^"^*^, we have 

< min{Cnl/2(log,^)(l + \^/k- ^^z\)-^ ,n^''^-'}- 
\Rtk\<Cn'"-'\ 
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provided \ Vk — E-/n\ < for some < 77 < 1/2, and 

\RTi\ < ni/2--/8A;-i/4ri/4|yfc _ 

ioY k and > n~^/^+'^. 

The proof relies on the well-known Plancherel-Rotach asymptotics for Hermite polynomials. These comprise 
three asymptotic expressions, corresponding to the the behavior of Ek in three regions defined relative to 
±^/2k, the order of magnitude of the largest zeros of Hk. In the case of Hermite polynomials, a classical 
reference is [15], Chapter 8. In [6], analogous formulas are derived for a general class of orthogonal 
polynomials. 

Theorem 5.1. Let Ek he the k-th Hermite function, then for any < /i < 1 and x in the "oscillatory 
region'': 

{x : Vfc(-1 + fi) < X < Vk{l - ^)}, 
we have the following asymptotic formula ((2.20) in see also (8.22.12) in flSf): 



2 , _i 7,1 x'^ 



-1/4 



cos 



Vkn I /9sc(2//^/fc) dy + - arcsin(x/\/fc) ) \ ^ + 0\\ 



+ sin(Vfc7r / psc(ty/VA;) dy — - arcsin(a;/"\/fc) I O [ — 
where psc denotes the semicircle density on the interval (—1, 1). 



(27) 



For x in the "transition region": 

{x : Vk{l - fi) < X < Vk{l + /i)}, 
we have the uniform asymptotics (JBjl, (2.21) and (2.23); see also (8.22.14) in JlSf ): 

1/4 



1/4 



Vk 



1/4 



1 



-1/4 



Vk 



-1/4 



A(x/Vfc)) Ai(/fc(:E/Vfc)) 



1 + -^] ' {fk{x/Vk)) '^^M'(^fk{x/Vk)))(l + o(r 



(28) 



when Vk{l — ^) < x < 1. If 1 < x < •\/fc(l — n), the same asymptotic formula holds with the overall sign 
reversed. Here Ai is the Airy function and Ai' its derivative and 

f-f (^))3/2 ^ i~kS^Ii Psc{y)dy l-fi<x<l. 
V -Iky)) j^y^TTiyTT^dy l<x<l + iJ,. 

Similar asymptotics hold in the transition region around the left edge: 

{x : Vk.{-1 - ^) < X < Vk{-1 + fi)}. 
Finally outside the transition region, for \x\ > Vk, the Hermite functions have Gaussian decay. 
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From the alternate form of the asymptotics m the transition region found in [TB], (see 8.22.14 there) 
combined with (|27)) it follows that: 

\Ek{x)\ < Ck-^/^\ 

for some C independent of k. 

We shall also use the following asymptotics for the location of the eigenvalues of the zero temperature case, 
see [6], Theorem 2.29 : 

Theorem 5.2. Let \k,n be the k-th zero of the n-th Hermite polynomial, where ko < k < n — ka for some 
ko, then we have: 



where is the inverse function to x ^-^ Psciv) dy. 



C 

- nH^{l-^)r/^ ^^^^ 



In this section and the next, we shall use heavily the fact that the eigenvector matrix for A„ oo is given by 
columns of the form 

1 

-(-E'„_i(Am), i?„_2(Am), ■ • • , Ei{X,n), Eo{Xm)) 



where A,„ is the corresponding eigenvalue and Ef^ is the fc-th Hermite function. The proof of this fact uses 
the three-term recurrence relation for Hermite polynomials. It can be found in or [S]. We shall also 
adopt the convention that the eigenvalues Xm are scaled to be of order ^/n and Am to be the scaled version 
that are of order 1. 

Proof of Proposition \4-S\ We will first establish the bounds for the diagonal entries Rkk ■ 

\Rkk\ < C7ii/2(log7i) • (1 + \Vk - V^?ftz\)-\ (30) 

First we write 

■m ™ 

where Um is the eigenvector corresponding to the eigenvalue Am, since Um{k) — ^ the above 

expression is 

^ _ Ek{XmY 1 

^ nEn-l{Xm)^ Xm - Z 

For the denominator, we have the following asymptotics [TB], section 2.6.4: 

nEn-liX,nf = nl/2(p,,(Am) + 0(71-^)), (31) 
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provided |Am| < 1 — 9(5/10. Since —1 + S<^z<l — S,we have: 

\Rkk\<Csn-'^' ^ ^^^ + Cs J2 

|A„.|<l-9<5/10 ' ^' |A,„|>l-95/10 

<C,n-i/2 V + (32) 

By the Plancherel-Rotach asymptotics (Theorem IS.ip . we have 

£;fc(a;) < Cmin|fc-i/4 j^l- ^ ,fc"'/''|. (33) 
Let E — , without loss of generahty, we shall assume that E > 0, and that \E — yfkjn\ = h/ ^/n^ i.e. 

^^^;| = h, 

with < /i = ft.(i?) < We shall consider the case where {) < E < yjk/n. The case where E lies in the 
exponentially decaying region of Ek{x) is simpler and follows by a similar analysis. 

First, note that the estimate ([50)) is efficient only when h > I. Consider the partition of the region 
(E, y^k/n) into the intervals 

Up = (v/fc?^ - pk-^^^n-^/^, ^fkf^ - (p - \)k-^l'^n-^l'^), \<v< hk^^"^. 

Note that each Up is an interval of length k^^l^-nT^I"^ . By the local semicircle law for the zero temperature 
case, the number of eigenvalues in an interval IJp is bounded by 

provided \Jp is at distance greater than c > from 1. Moreover, we have 

Eu{xf < Ck-^'^'p-^''', X e Up, 

and 

I A™ -z\> \hn-^'^ - pk-^"'n-^'^\ + n-^'^+\ 

So the contribution to the sum in equation ((5^ due to eigenvalues Am lying in the region {E, ^Jk/n) is 
bounded up to a constant by 

fc-l/6„-l/2 ^ fc-l/6p-l/2(|;,^-l/2_p^-l/6„-l/2|^^-l/2+.)-l 

Simplifying the above sum, and letting h — k"' for some < a < 1/2, we have the bound 

J2 p-^/^{\l-pk-^/^-°'\+n''k-'')-\ 

p=l 
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The sum over p is bounded by C(log and so the above quantity is bounded by Clogn • 
^_i/4-a/2 < j^-i^ gjj^pg < _Q,/2. 

Turning to the contribution to the sum in p2p from eigenvalues A„i in the region (— we have that 

Eu{xf < < h-\ 

for h < fe^/^. Combined with the fact that ^ < Cnlogn (see Lemma [6.41 below') . this establishes 

the estimate ([50)) . The alternative bound |i?fcfc| < n^/^~'^ follows from Cauchy-Schwarz inequality, noting 
that the eigenvectors have norm 1. 

Finally, to prove the bound for Rkk when h < , < r/ < 1/2, note that in this case 

k ^ En + 0{n^/^+'^). 

The result then follows from ([33| and ([32]) by a dyadic decomposition as previously. As noted above, the 
implicit constants depend on rj, but we only used this bound for the fixed choice rj = 1/4. 

We will continue by proving the bound for Rjk- To simplify notation, we shall assume that E = — 0, 
the case for nonzero E follows similarly. Unlike in the case of the diagonal entries Rkk, we will exploit 
cancellation. 



We write Rki as 



{k)u,n{l) 

Rki = 2^ — 



T • "(ft 

m— 1 



Applying summation by parts, we rewrite the above as 



n-l . s n 

E ^ \ ]Y.^t{k)ut{l), (34) 

m— 1 ^ ^ 



t—7n 



Note that: 

EkiXt)Ei{Xt) 



t=m t=mn{En-l{Xt)) 



From here on, we shall assume without loss of generality that k > I. We will concentrate on the contribution 
to ([55]) from eigenvalues |At| < (1 + The sum over A* lying in the region where Ei{-) decays 

exponentially is simpler to handle, and uses similar techniques. The following is a key lemma: 

Lemma 5.3. Given 1 < I < k < n, let m be such that: 

Xm e {-V1{1- fi),Vl{l- fi)), 

for some fixed < jj, < 1 independent of k,l,n. Then, for any m such that 

XrH e i-Vl{i- pi),Vl{i- 
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we have: 

Ek{\t)Ei{\t 



t—m iL ± \ J 



< cfc-^/^r i/4(yfc - Vi)-\ (36) 



\Xt\<il-^^)Vl<{l-^i)Vk, 



Proof of lemma \5.S[ Since 
equation (j27p implies: 

EkiXt) = ^k-^^^(^l-j-^ cos (^VkirJ^ p,,(y/Vl) dy + ^ arcsin(At/\/fc) + j ^ ^ 
+ sin(VkTT J psc{y/Vk) dy + ^a,rcsm{Xt/ Vk) + ip,n,kj O (^^^ , 



where 



<^m,fc = / psciy/Vk)dy 

= klT J Psc{s) ds. 

A similar expression holds for Ei{\t). For m fixed, we will expand Ek{Xt)Ei{Xt) and (£'„_i(Ai))^ around 
Am- By ([29| . we have 

A* = (^C (^^^^ + ^^^^^ + ^ arcsinC((TO + (i - m))/n)^^ + 0^(^-2). 
By Taylor expansion around Am, we have 

At = Am+ai^^ + 0^(n-3/2), (37) 

for < t — m < ri. The constant 

,,/6m-3 1 . , ^^ / 1 C'W") \ 

ai = C —:. K 7: — arcsinC(m/n) 1 + v / / 1 



is independent of t. \ai\ is bounded above and below, uniformly for m such that A„i G (^(1 ^ (1 ^ 

fi)^/n). Combining the above and using the product to sum formula cos a cos b — (1/2) •cos(a+&)4-cos(a— fo), 
we have the estimate 

EkiXt)Ei{Xt) = Ck-^^H~^/^cosiVk7: [ 'p.dy / Vk) dy - VIt: [ ' psciy/Vl)dy 



Am "'A 



+ i arcsin(At/Vfc) - i arcsin(Af /\/I) + ipm^k - ^mj) ^1 + O 
+ (cos(-) term with sum of arguments) . 
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The final term denotes an expression formally similar to the initial part of the right side of the equation, 
except that in the argument of the cosine, all — signs are replaced by + . We will denote ipm = ^m,k — ^m,i] 
this term is independent of t. Taylor expansion of the integrals and arcsin functions in the expression above 
yields: 

Ck~^/H'^^/^(^cos (n{^/kpsc{Xm/Vk) - ^lpsc{\n/^l) + ^(arcsin(A™/%/A:) - arcsin(A„/v^))) • (At - A^) 

+ 0((%/fc - ^^l)/^/k) ■ (At - \r,^f + ) + O 

The constants implicit in the O(-) terms are bounded uniformly in m such that Am G ((— (1— 
Using the expansion (I37p of At around A„i, we are lead to the following approximate expression for the 
product Ek{\t)Ei{\t): 

C;,-i/4pi/4 (^ai{^ - Vl) ■ 

+ Ck-^'H-^'^ ■ (oiiVk - Vl)/Vk) ■ + 0{{{Vk - Vi)/Vk) ■ n-\t - m)2)) , (38) 

where ai is the constant independent of t in p7l) . We have suppressed lower order terms in the final 
displayed equation. In particular, we have omitted the second cosine term referred to above, involving 
sums instead of differences in the argument. When t varies, this term oscillates faster, and the remainder 
of the argument will make it clear that it is bounded by the contribution from the term in the last equation. 




We wish to use the approximations we have introduced to estimate the sum in ()5.3|) between A™ and A^. 
Both eigenvalues are assumed to lie in the range ((— 1 + fJ-)Vl, (1 — fJ')Vl), so that \Xm — Am| < 2Vl. Thus 
there are 0{Vnl) indices t such that Am < At < A^. To take advantage of the oscillation of the cosine 
term in our approximation, we perform the sum "period by period" . 



For fixed m, we sum Ek(Xt)Ei{Xt) /n{En^i{Xt))'^ over indices t in the range 

m <t < m -\ j= -j=- = m + A, 

ai{\k — V 

corresponding to a full period of the cosine term in (p8| : 

^ Ek[Xt)El(Xt) , ^ , s, 

-— — — = I(m + II m . 

n[En^i{Xt)y 

t— m— \ 

Here, I(m) represents the contribution to the sum from the first (oscillatory) term in psp . The sum II(m), 
the contribution from the error terms, is readily estimated; by (I38p and (j3ip . we have for any m: 

|II| < Cfc-l/4;-l/4(^-3/2^-l/2 + (Vfc _ Vir^k''/^). (39) 



To estimate the sum I(to), we first consider the case A > 100. After using a trigonometric identity, we are 
faced with sums of the form: 

c(m) J2 '^"^ T^) ' 

k=0 ^ ^ 
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where c{m) depends on the phase if{m), and |c(m)| < 1. The sum is bounded uniformly in A. Moreover, 
the sum remains bounded when the upper hmit of the sum is varied by less than 5: 



LAJ±r 



k=0 



< C, r < 5. 



By choosing < r < 5 appropriately the expression in (|40[) can be made either positive or negative. There 
are 0{Vl{Vk — Vl)) periods separating Am and Am- We apply the above approximation repeatedly, using 
the bound (|39p . and adjusting the upper limit of the sum over each period so as to make the signs of I 
alternate. The estimate ((36|) thus follows in the case A > 100. In case A < 100, we must modify the 
previous proof somewhat. The error terms in (j38l) can be treated as before, but when dealing with the 
oscillatory term, we are no longer assured that a small change in the upper limit will change the sign of 
the trigonometric sum in (j40p . The phase term tpm must be taken into account. When m' > m, we have: 

Vm' =Vm. + VA-tt / Psciy/Vk) dy - VItt / Psciy/Vl) dy. 

= + 2{Vk - Vl)(Am' - Am) + O ^^^=-^(Am' " Am)^^ . 

Using this expansion when m and m' define the limits of two consecutive periods, we can ensure that the 
oscillatory term in (j38p . when summed over different periods, remains bounded by Cfc-iZ-^r 1/4^^-1/2. The 

□ 



approximation of the ipm' introduces an error that is controlled similarly to the terms in ((39 



In the transition region, the behavior is similar: 

Lemma 5.4. Given 1 < I < k < n, let m be such that \„i G (— %/I(l — fJ,),Vl{l — /i)) for some fixed /i 
independent of k, I, n, then there exists rh such that Am G (— VI(1 — /i), — j-i)) 

< Ck-^'H-^/^n'^iVk - Viy^ (41) 
for any rj > 0. 

Remark 2. We will only prove the lemma for the right edge. The proof of the statement for the left edge 
at —1 is identical. 



t—rn ^ \ / / 



E 



Proof of lemma \5^\ The proof uses ideas similar to those in the proof of (j5.3|) : we replace Ei and Ek by 
appropriate asymptotic expressions and thus reduce the problem to estimating some trigonometric sums 
with slowly varying frequencies and phases. We will omit some details and give the orders of magnitude 
of the approximations and error terms involved. 

We use asymptotics for the Airy function and its derivative (see pi, Ch. 10, section 4): 

Ai{x) = -^^^e-i-'^^(l + 0{x-^/^)) (42) 
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and 

We first establish an estimate for tlie magnitude of Ei in the transition region. If A = l^/'^ — /-i/2+4a fQj. 
some a, then have 



-1/4 



Using the asymptotics for the Hermite functions in the transition region pS)) and the asymptotics for the 
Airy function for x close to A, we find: 

Ei{x) = Cr" sin I — psciv) dyl\ 

= sin (^(p..(A/^/^)) • ^ + 0(p;e(A/VI)) • ' ) . 



We first consider the situation where y/k — \/l > jiy/l. Fix an eigenvalue A„i with 

Near Am, we have: 

Ek{Xt)Ei{Xt) = k-^/H-^/^ cos I Vkn f ' Psc{y/Vk) + ^ arcsin(At/%/I) 



X 1 + 



1/4 



1 - 



hi 



-1/4 



(44) 



{f{Xt/v^))-'/'A^{f{Xt/v^)). 



Here again we omit lower order terms, as they are bounded by the expression above. Write 

A„^Vl-rV2+4«, 

for some < a < -j and ^-1/2+4" ^ pl^/^. By a previous remark, the magnitude of Ei around A,„ is 
Cl^°' . For At near Am, the approximate expression for Ei found above becomes 

El{Xt) = cos (^P"{Xt - Am) + 0{l'/^-^"{Xt - Am)') + V-t.m) , 

where ipt,m denotes a phase term independent of t. From equations (jH]), (|29p and ([3T|) . it follows that 
there exists rh of order m + 0{y/n{k^/'^ — P°')~^) such that the following estimate holds: 



E 



EkiXt)Ei{Xt) 



n{En-iiXt)y 
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Consider a partition of the interval (Vl — I i/^+^a^ _ ; i/e^ [^-^q subintervals of the form 

(Vi - r VI - r i/2+4"-pr()^ p = 1, 2, . . . . 

Here 77 > is some parameter to be determined. There are 0(77^^) such intervals. The length of each 
interval is ^-i/2+4a ^^^^^ ^-j^g sum over Am belonging to such an interval is bounded, up to a constant 
factor, by 

^-l/4^-a^-l/2+4Q^l/2-2Q^^l/2 _ fay2 

This last expression is Ck^^^'^l^^^'^{Vk — Vl)"-'-. In summary, for any A^, such that A™ < (1 — l^)Vl, there 
exists some Am near \/I — l^^/^ such that 

The contribution to the sum from eigenvalues A* in the region between \/I — /^^/^ and Vl can be trivially 
bounded by Cr i/i2/-i/6)!j-i/4^-i/2 

using the amplitude bound for Ei in the region and the oscillation 
from Ek ■ This establishes the case where — a/I > fj,Vl. 

Next, we consider the case where 

for some /3 such that < fil^^^. In this case, the sum involves contributions from products Ek{Xm)Ei{Xm) 
such that Am lies in the transition region for both Ei and Ek- Suppose 

Am = x/I-rl/2+4a, 

for some a such that \Ei{Xm)\ — 1^°"- Since k > l,we have the following estimate (writing m = |^-)-^-i/2+4a 
to simplify notations): 

Ekix) = Ck-'^/^m'^/^sm (k^/'^m^/^{x - A) + 0{k^/'^m-^/^) ■ {x - . (45) 

A similar expression holds for Ei. Thus, for Af near Am, the product Ek{Xt)Ei{\t) can be replaced by a 
sum of trigonometric functions oscillating at frequency fc^/'*m^/^ — P", at the expense of introducing an 
error of order 

0[k^'^m-^/^ -l^'^-^")-(x-\f. 

Using the error estimates ((29)) and (1451) and summing over the eigenvalues in a region of size 0(Z^i/2+4a+'7) 
(containing 0(/~^/^+''"+''(fc^/'*m^/^ — 1'^°')) periods) will introduce an error of 

1,1/4^-1/2 71/2-20! 
rfr-l/8rT7-l/4;-l/2+3a+r, . ^ ~ ' 

(fcl/4^1/2 _ ;2a)2 • 

We consider two cases: (i) /3 > — ^ + 4q and {ii) /3 < — + 4a. For the first case, the above expression 
simplifies to (recall the definition of m): 

(7;t-l/8-l/4;-3/3/4-l/2+3a+77 . |-^l/4^-/3/2 _ ^l/2-2a-) < ^^-1/4^-1/4^-/3+,)^ 



23 



In case (n), we have the bound: 

(7;(.-l/4^-l/8^1/8^1/4-4a^,,^ 

Since a > ^ + j, the above expression is bounded by Ck^'^/^l^^/H'^T' . Finally, since = (\/fc — 
proceeding as in the case of — \/l > /Lf\/I, we arrive at the desired estimate. □ 



The contribution from eigenvalues Am in the region (Am > Vl + l^^^^^'^) where Ei decays exponentially is 
easier to handle than the cases treated above. The case Vl < \Xm\ < Vl + 1^^/^+*^ can be estimated in a 
similar manner to the case (1 — < \Xm\ < Vl- Combining the two lemmas above and first summing 

over Am in the two transition regions, and then the oscillatory region, we have established the following 
proposition: 

Proposition 5.5. For any m, we have the hound 



n{E 

for any rj > and some constant C . 



We substitute the estimate in the previous proposition into equation p4p to obtain the following bound 
for Rki (recall that we are assuming 3?z = 0): 



\Rki\ < Cnn-'/^k-'/^iVk - Vl)-' J2 



1 



Am ^ Z Xrn+l ^ Z 

To estimate the above, perform a dyadic decomposition around !Rz. Consider the eigenvalues Am such that 

2P„~l/2+. < |Am| <2f+l71-l/2+^ 

Their number is 0(2^'n~^/^+'^) and the factor | ^ ~ >, I is 0{n~^'^2~^P). Therefore, the contribution 
to the sum from such eigenvalues is bounded by Cn^/^^'^2~P. Summing over p, we find: 

\Rki\ < Cn^n''^-H-''^k-'/\Vk-Vl)-\ 

Taking, say, 77 = e/10, we have established Proposition CI 

6 From 1/2 to 1: Inductive arguments 



In this section, we improve the semicircle law from the level of n 1/^+'^ to the optimal level of n ^'^'^ by 
an inductive argument. 
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Theorem 6.1. Let s{z) be the Stieltjes transform of the measure induced by the eigenvalues of the nor- 
malized Gaussian (3-ensemble, -^An^fs- Let Ssc{z) be the Stieltjes transform of the semicircle law. Then 
with overwhelming probability: 

sup \s{z) - s,c{z)\ ^ o{l) (46) 

where for e > and < S < 1, the domain D^^s is defined as 

:= {z : 3z > n-^+\ -I + 5 <^z < I - 5}. 



Before we start we will need two facts about the tridiagonal models and Stieltjes transforms. The first can 
be found in [7 : 

Proposition 6.2. The tridiagonal model Ap can be diagonalized as A = QAQ* , such that the first row of 
Q is independent of A and consists of independent entries of xp- distribution normalized to unit norm. 



The second lemma establishes the link between control of Stieltjes transform and control of distribution of 
eigenvalues and can be found in 17 : 

Lemma 6.3. Let 1/10 > rj > 1/n and L,e,T > 0. Suppose that one has the bound 

|s„(z) - s{z)\ < T 

with overwhelming probability for all z with |5Rz| < L and 3z > rj. Then for any interval I in [—L + e, L — e] 
with \I\ > max(2?7, ^ log one has 

Nj-n j Psc{y)dy <e Tn\I\ 
with overwhelming probability, where Nj denotes the number of eigenvalues in /. 



Let us introduce some terminology: we will say the semicircle law holds at level a if, for z G Di^a.s, and 
any c > 

\sfi{z) - Ssc{z)\ < C 

for sufficiently large n, with overwhelming probability. We first prove the following simple lemma: 

Lemma 6.4. Suppose the semicircle law holds at level Qz > n'^ (or equivalently by lemma [U73\ the Stietljes 
transform is close to that of the semicircle law at the level n"" ) for some — 1 < a < 0, then we have 

iE^<^(^^)-^""logn (47) 

for any z such that < 5z < 
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Proof. We will in fact show the following inequality: 



The lemma then follows from the trivial bound: 



|Aj - z 



holding for all j. Suppose the semicircle law holds at level n ^ < 77 < 1/10, that is: 

\s{z) - Ssc{z)\ = 0(1) 

for any z such that 5z = 77. First notice that we have the following trivial bound for any i 



3 

where 3z = 77. It remains to show that 



^ \X,-w\ - Qw^\X,-z\-^ 



\Xj ~ z 



for such z. Let E ~ and consider a dyadic partition around E: 

Up := {j : r-^Ti < \X, -E\< 2^ 

for p = 1, . . . , log r/~^ 

Uo:^{j :\X,~E\<T^}; 
Uoo ■.^{j:\X,~E\ > 1}; 
By lemma 1131 < 2Pnr] and over the set Up, — 2^~''?7^^, thus 

Summing over n, we obtain the desired bound after noting that log 7/"^ < logrT,. 
The proof of the theorem relies on the next two propositions. 

Proposition 6.5. Suppose the semicircle law holds at level a for some — 1 < a < 0, then 

|<(z)-S,e(z)|=o(l) 

for z such that 3z > 7j(i-i)/2+<5 Jqj. ^j^y J > Q with overwhelming probability. 
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Proposition 6.6. Suppose \Rii{z) — Ssc{z)\ = o(l) for z such that Q^z > for some —1 < a < 0, then 
we have an improved semicircle law, i.e. 



n ^-^ A," — z 

3 ■' 



o(l) 



for z such that 3z > n'^"' ^^^^^ for any S > with overwhelming probability. 



Proof of Theorem \6.1[ By Proposition l4.1[ we have a local semicircle law at the level of —1/2 + e. Applying 
Propositions 16.51 and 16.61 repeatedly will yield the desired conclusion. □ 



What remains is the proof of the two propositions above: 

Proof of Proposition 1 6. 51 By Schur's complement, we have the following relationship 

1 



^"11 „/3"12-' 

where an is the (l,l)-entry of our (symmetric) tridiagonal matrix A, normally distributed with mean 
and variance 1 and 61.2 is the (l,2)-entry of A, which follows a ^-distribution with (n — degrees of 
freedom. R is the resolvent of the matrix obtained by removing the first row and first column of A. We 
will denote the latter matrix by A. 

Firstly, notice that 011,612 and R are independent, and that -7i=|aii| < n^^/^^"^ with overwhelming 
probability for any c > 0. Thus, we have 

^ aii-z-^bl^Rii = -z-^Ebbl^Rii + ^{bl2-mi2)Rii+Oin~^/^+^). (49) 



y/np np np np 

For the second term on the right hand side, we have £(,612 = {n — l)/3, so 

—Ebbi2Rii — -Rii- 

np n 

Let ((71, . . . Qn-i) be the first row of the eigenvectors for A, Xj be the eigenvalues and write 



A," — z 

3 3 



Similarly, 



if 



A, — z 
3 ■' 
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We would like to compare Rn and Rn, to this end we write 

Rii Rii = ^qy2 — 



— —--Z^J-— (50) 



By proposition 16. 2[ q and A^ are independent and so are q and Aj, 

o2 o2 _ i 

y ^^-E.y ^^ = y — ^. 

^ A,- - z A, - z ^ Ao - z 

3 3 3 

Since q and Aj arc independent, we can condition on Aj and apply McDiarmid's inequality (see the ap- 
pendix) to the sum to conclude that the term on the right hand side is bounded with overwhelming 
probability by 

1/2 



,1-c |\ , _ 



By Lemma [6.41 this last quantity is in turn bounded by C(3z)^^rt^^/^+''/^+'^ logn for some constant C 
and any c > 0. The fifth and sixth term of equation ()50p are similarly bounded. For the middle two terms, 
we used the interlacing property of the eigenvalues of a matrix and its minor, to obtain that the term is 
bounded by — , where ?7 = Qz. 



The last thing we have to control is the term -^{h\2 ~ ^bb\2)Rii in equation Since h^-s the x- 

distributions with degree {n — ;^(&i2 — ^^6^2) is bounded by with overwhelming probability 

for any c > 0. By equation, (gg]) andjl^l -Rii < C(5z)"in°+'= for any c > 0. 



Putting everything together, we have the following relationship: 

i?ii(z) = + 0((3z)-ln-l/2+'^/2+'= + {wq)-^ + (3z)-ln-l/2+a+c)^ 

-z - Rii 

The condition of the proposition guarantees that the error is o(l), and by the standard argument of 
inspecting the functional equation of the Stieltjes transform of the semicircle law, the proposition is estab- 
lished. □ 



Proof of Proposition \6. 
establish that 



The proof is similar to that of Proposition 16.51 By the assumption, it suffices to 



n A, 



-<(z) 



0(1) 
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with overwhelming probabihty. The difference inside the absolute value sign is ^ 
concentration of measure argument, and lemma the statement is established. 




— . So, again by a 



□ 
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8 Appendix 



Here we prove a slight modification of McDiarmid's inequality used in our inductive argument. 

Proposition 8.1. (A modification of McDiarmid's inequality.) LetXi,...,Xn be independent subgaus- 
sian random variables and suppose F is a function of n variables such that there exists an event with 
overwhelming probabibility that if xi, . . . , Xn, £i € ^, then 



for all 1 < i < n, and outside of Q, F is bounded by a polynomial of n. Then for any A > 0, one has 



for some absolute constant C,c > 0, and a — J2i=i 'hi- 
proof. By symmetry, it suffices to show 

F{F{X) - EF{X) > Act) < CexpicX^) 
Let i > be a parameter to be chosen later. Consider the exponential moment 

E{exp{tF{X))\Xi, . . . ,X„_i, {X,} e n). 
Writing Y = F{X) - E{F{X)\Xi, . . . , X„_i, {X;} e fl), we can rewrite the above as: 



By the condition of the theorem, tY fluctuates only by at most te„ and has mean 0. By Hoeffding's lemma, 



\F{xu-.-, 



P(|F(X) - E(F(X))| > Act) < Ccxp(-cA2) 



E(exp(tF(r))|Xi, . . . ,X„_i,X, G n) exp(tE(F(X))|Xi, . . . ,X„_i, {X,} e ft). 



we have 



E(exp(iF(r))|Xi, . . . ,X„_i, {X,} en)< exp{0{t^cl)). 
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Integrating out the conditioning, we have the bound 

E{exp{tF{X)) < eMO{t^cl))E{t{EF{X)\X,,. . . ,X„_i, {XJ e fl)). 

Now the latter expectation is a function of the first (n — 1) variables and obeys the same hypothesis, so we 
can iterate and obtain the bound: 

n 

exp(^O(t2c?))E(t(EF(X)|{X0 e n)) 

i=l 

and by the overwhelming probability condition and the bound outside of the set Cl, we have 

¥{F{X) - EF{X) > Xa) < exp{0{fa^) - tXa) 
Optimizing in t gives the desired result. □ 
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